Quantum phase transitions of magnetic rotons 
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Due to weak spin-orbit coupling, the magnetic excitations of itinerant ferromagnet become mag- 
netic rotons, excitations with degenerate minima on a hypersphere at finite wavevector. Using 
self-consistent Hartree and renormalization group calculations, we study weak fluctuation-driven 
first-order quantum phase transitions, a quantum tricritical point controlled by anisotropy and non- 
Fermi liquid behavior associated with the large phase volume of magnetic rotons. We propose that 
magnetic rotons are essential for the description of the anomalous high-pressure behavior of the 
itinerant helical ferromagnet MnSi. 
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The theory of classical, second-order phase transitions 
has been extended to quantum phase transitions^, es- 
pecially so in itinerant electron systems 0, 01 • Critical 
non-Fermi electron-liquid properties and new phases can 
emerge near the quantum critical point, and unconven- 
tional behavior of quantum criticality attracts a lot of 
current interest 0, |^ IE S 0- Among different scenarios 
for unconventional behavior near the phase transition, 
a generic scenario associated with a large phase space 
of fluctuational modes has been little explored^, even 
though it is well-known for classical phase transitions [lo|. 
A classical phase transition becomes fluctuation-driven 
first-order if soft fluctuational modes have a minimum 
at a nonzero wave vector |q| = go on a hypersphere in 
il'-dimensional space |loj. In the presence of spin-orbit in- 
teraction such soft modes with large phase space appear 
naturally (see below) in the proximity of a continuous fer- 
romagnetic transition and can be called magnetic rotons 
(in analogy with roton excitations of superfluid He^). In 
this letter we present the theory of quantum phase tran- 
sitions of magnetic rotons. 




FIG. 1: The lowest paramagnon mode of a helical magnet 
becomes, for weak anisotropy, a magnetic roton with minima 
for all |q| = go > 0. The large phase space of magnetic roton 
fluctuations is responsible for the unconventional behavior. 



We show that the large phase volume of magnetic 
fluctuations with minima at |q| = q^, independent on 
the direction of q, is responsible not only for the T^'- 
order nature of the transition and strong non-Fermi liq- 
uid damping of itinerant electrons, but also for a new 
tricritical point controlled by anisotropy. We combine 
the theory of quantum critical fluctuations, as developed 
in Ref.0,l3|j with the theory of weak "crystallization" via 
fluctuation-driven first-order transitions, as developed in 
Refs.|E 0, ^3- In addition, the resistivity of itinerant 
helical ferromagnet and the electron self-energy are cal- 
culated due to the scattering of electrons by roton fiuc- 
tuations as well as the strong tendency towards glassy 
behavior. Finally, we argue that magnetic rotons form 
a starting model for the description of the anomalous 
high pressure behavior of the itinerant helical ferromag- 
net MnSiQ. 

We introduce a collective magnetization M (z, r) along 
the lines of Refs.0, 0, 0|, valid for a crystal without 
inversion symmetry, like MnSi, and obtain the effective 
action S [M] = S'dyn [M] -t- S'stat [M] which consists of a 
dynamic term (see below) as well as the static term 

5,tat [M] = j d^^dr (^^ |M|' + J |M|^^ 

+ J d^xdr ^1 ^ (VAf,)^ + 7M • (V X M)^ . (1) 

Here, i = {x,y,z) and S determines the proximity 
to the transition point. The last term is due to the 
Dzyaloshinskii-Moriya interaction, caused by spin-orbit 
interaction in non-centrosymmetric crystals. The con- 
stant 7 <^ a/a is small due to the relativistic origin 
of the spin-orbit coupling. We assume analyticity of 
the Landau expansion, Eq|Tl and a,u > 0, while other 
scenarios were suggested jlJQ, In Eq^ we have 

not included anisotropic terms like a' J2i {di^if and 
u' ^ j Mf allowed by crystal symmetry, but smaller by 
[q^aY <C 1. Quantitative estimates |0| demonstrate that 
the Dzyaloshinskii-Moriya interaction is indeed the dom- 
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inant effect due to spin orbit interaction in MnSi. 

Spin-orbit coupling (7 7^ 0) splits the otherwise de- 
generate paramagnons into longitudinal and transverse 
modes (see the Fig.l). The energy of the lowest mode, 
E {q) = Aq -I- (|q| — qo) /2, has minima for all |q| = qo, 
i.e. is a magnetic roton. The mode corresponds for 
7 > to a right handed helical mode. The bare gap 
Aq — S — aqQ/2, which might be tuned by varying 
pressure or chemical composition, will be renormalized, 
Aq A, due to interactions. If A is much smaller than 
^qQ, the energy difference between the roton and the 
other modes, those modes can be omitted for ksT < ^qg. 

The dynamic part of the action depends on the cou- 
pling between magnetic rotons and particle-hole exci- 
tations. If the roton gap A is inside the particle- 
hole "stripe" in (oj, q)-space, the damping is linear 
in frequency 2]: S'dyn = /^^|M|^. Here = 

T J2n I ^^^'^ ^ ~ ^ with fermion-roton couphng 
constant, g. Since in general, it is complicated (espe- 
cially in ordered helical states) to describe the particle- 
hole excitations '17^ and thus the spectrum and damping 
of rotons, we introduce, for simplicity, the variable dy- 

namic exponent z via \lu\ — s- Ep ( ■g^j . z = 2 cor- 
responds to the over-damped case, while 2: = 1 to the 
undamped case. As we will show, distinct universality 
classes occur for different z-values, making it natural to 
allow z to vary. In fact, e — 3 ~ z plays a role simi- 
lar to the D — A expansion parameter in classical critical 
phenomena. 

Considering solely fluctuations of the lowest mode we 
obtain the effective low energy action 



S = - I Xq "Vq'P-q + 



'Pqi^q2^q3^q4 (2) 
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with qi = —qi — q2 — q3 and 



X 



-i(ro) = |c.p/^ 



(|q|-i) 



(3) 



expressed in terms of dimensionless variables: q — > q/^o 
for momenta, tq = Ao/(q!(7q/2) for the roton gap, to 
-^{2Ep /Taq^y/'^ for frequencies and A 
for the coupling constant. Notice that for D ^ i the 
coupling constant A <C 1 is small for z > 1, making 
a weak coupling expansion much better defined than 
in the classical limit[10j. The momentum dependence. 



A, 



qi,-qi,q2,-q2 = I + (qi • q2)^) results from the pro- 
jection onto the roton mode. It is useful to introduce 
two coupling constants: A averaged over the angle be- 
tween qi and q2 (i.e. for generic angles) and A|| where 
all wave vectors are parallel For a momentum inde- 
pendent initial interaction, A|| = |A = A. A weak mo- 
mentum dependence of u can change A/A|| further. Many 
anomalies related to rotons stem from the fact that, at 
low energies, the integrals over momenta are essentially 



one-dimensional and that A|| renormalizes qualitatively 
different from A. 

The nature of the phase transition can be understood 
from an analysis of the roton gap r and the interac- 
tion vertex \{qi). A self-consistent equation for the gap 
within Hartree approximation, valid for small A, can be 
obtained from the Dyson equation for the field propaga- 
tor: r = r{) -\- XTiH{r), with singular Hartree self-energy. 



(r) ^ Xq{r)^z 



A' 



27r2 (z - 1) 



(4) 



proportional to the local fluctuations of the magnetiza- 
tion, and physical roton gap, r. A is the upper momen- 
tum cutoff. For < z < 1 the fluctuations of the lo- 
cal magnetization diverge as r — > 0, like for the classi- 
cal (T 7^ OA') fluctuation-driven l''*-order transition ex- 
plored by Brazovskii lOj. In contrast, T^h [r 0) stays 
finite for z > 1, relevant for itinerant magnets, and no 
disordered state is allowed for rg < — AE^f (0). In the 
ordered state, the gap of magnetic rotons 



r = ro + AEH(r)+6A|mop 



(5) 



is obtained by adding the Hartree energy due to an helical 
order parameter with amplitude mo. toq is determined 
by the equation of state 



mo?" = 3A||Kmo^ 



(6) 



with K = 2A/A|| — 1. For ro smaller than the value 



^spino a _ — ASh (0) + CzX (Ak) and z < 3 three ex- 
trema of the energy exist if k > 0. is of order unity. 
Decreasing ro, the lowest energy solution jumps from 
7710 = in the disordered phase to 7710 7^ in the ordered 
phase. The roton mass r discontinuously jumps to a 
larger value in the ordered state. The regime of metasta- 
bility, ro^spinodai + ASff (0), vanishes as k — > 0. For k < 0, 
777,0 = is the only solution of the equation of state, but 
solutions exist for several ordering directions of qo (see 
below). For 1 < z < 3, the phase transition remains T''- 
order for k > 0, despite the fact that the Hartree energy 
does not diverge as r — > 0. It is important to notice that 
for r ^ the r-dependent part of the self-energy remains 
the dominant one: r(^~^)/^ » ro -l- AI]//(0) . 

A direct way to demonstrate the T^'-order nature of 
phase transition is to analyze the renormalized interac- 
tion vertex Ajj. The dominant renormalization of this 
interaction comes from the polarization "bubble" dia- 
grams. The polarization "bubble" is defined as n(r) — 

Xq {rf and behaves as H(r) = ~ r(^-3)/2. It 

is necessary to account for a constructive interference be- 
tween two channels of zero total momentum^CJj . Still, the 
perturbation series is dominated by terms due to H(r) 
(for a proof see below). The summation of the leading 

terms gives: A'| = -^llT^j^^^^- If n(r 0) diverges. 
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FIG. 2: (a) Flow diagram for the two coupling constants \a 
and Ai, in units of A*.(b) Resulting Aa-ro phase diagram for 
fixed At with quantum tricritical point at \a = 2Ai, separating 
the regimes with second and first order phase transition. 



AjJ changes sign for k > 0, implying a I'^'-order tran- 
sition. The renormahzed vertex A = — - , , with 

1+An(r) 

generic angle between momenta behaves differently, and 
does not change sign. The case 2: > 3, corresponds to 
an ordinary 2"'^-order phase transition above the upper 
critical dimension. In contrast, z < \ corresponds to the 
dimensionality below the lower critical dimension for sec- 
ond order transitions. Therefore we have to distinguish 
between two classes of fluctuation-driven 1-order transi- 
tions characterized by diverging (z < 1 "classical like") 
and non-diverging (1 < z < 3 "quantum") local fluctua- 
tions. 

In the two-loop approximation, an important differ- 
ence between an ordinary c^^-theory and the roton field 
theory becomes evident. Diagrams not taken into ac- 
count in the above renormahzed vertices. Ay and A , are 

small by a factor r'^"-'^'/^ which comes from the angle in- 
tegration if the total momentum of the propagators inside 
the "bubble" is not zero. This implies in case of A = A||, 
where the I'^'-order transition occurs if 1 — An(r) ~ 0, or 
7-(3-2)/2 ^ WmX self-consistency of the Hartree solution 
is guaranteed by the condition <^ 1, much 

less stringent than for classical systems 



qualitatively different. We obtain the flow equations: 
^ = £A-AV(Te-',r(0) 



dA| 

dl 



eA|| -2A /(Te-',r(0) 



(7) 



In the case A|| = 2A (k = 0) the system is right between 
a regime with a fluctuation induced first order transition 
and a regime where the equation of state gives mo = 0, 
suggesting the existence of a new critical end point or a 
tricritical point. Insight into the roton model can then 
be gained by using a renormalization group approach, 
following Ref. 11, 18]. As suggested by our earlier calcu- 
lation, the two coupling constants, A and A|| , renormalize occurs 



together with a corresponding flow equation for r(Z). 
Here / (T, r) = 3 Xq ('') with shell integration over 
momenta A/e' < |q — Qol < A. In agreement with 
the previous paragraph, our calculation is controlled if 
e = 3 — z is small. The zero temperature flow is shown 
in Fig. 2a and is characterized by an unstable fixed point 
A| = 2A* cx e. The fiow away from this fixed point is 
towards A|| ±00 and Aa — > A*. The fixed point can 

only be reached by tuning r and the ratio i.e. it is a 
tricritical point or a critical end point, see also Fig. 2b. 

The flow towards A|| —00 corresponds to a fluctu- 
ation induced first order transition[2(ij| . If A|| (l) changes 
sign at the same scale where scaling stops (r (l) = 1) 
we obtain the bare roton gap at the spinodal, ro^spinodai, 
identical to the result obtained within the Hartcc ap- 
proach. 

Close to the fixed point we can perform a scaling anal- 
ysis. The free energy of rotons to first order in e is 

F(f,A||,T,/i) =&-(l+^)F(f6l/^A||6^T6^/^6^'•) (8) 

where ~ 2 — je. feisa scaling parameter and r and 

A|| are the deviations of r and A|| from their fixed point 
values, respectively, /i is a field conjugate to the order 
parameter with = Scaling arguments yield for 

the correlation length ^ (r) ~ r^'^ , £, (^||) ~ -^n 

^ (T) ^ T~~ , depending on how one approaches the crit- 
ical point. Varying temperature at the fixed point gives 
C - T- for the singular contribution to the specific heat 
and Ti-i - for the NMR spin lattice relaxation 

rate. The discontinuity of the order parameter at the 

first order transition vanishes like m ~ A 
proaches the fixed point 
such that an additional w^p 

Finally, we discuss the flow A|| {I) —> 00, which occurs 
when the bare values obey A < 2A||. From the equation 
of state follows that mg = is the only allowed solution. 
However, for z > 1, the renormahzed roton gap vanishes 
for ro — > —AS// (0). If r becomes arbitrarily small there 
will be a crossover scale where the ignored anisotropy 
terms, no matter how small, dominate the low energy 
physics. States on the hypersphere are not degenerate 
anymore, and only isolated points of low energy excita- 
tions become relevant. The universality class is then the 
quantum version of the theory analyzed in Ref.|2ll| (see 
Ref . [liji ) . For Z? = 3 and z = 2 a mean field second or- 
der transition to a state with fixed direction of the helix 
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as one ap- 



AU these results require z > 2 
^ interaction is irrelevant. 
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A crucial difference of the two ordered states readied 
on either side of the tricritical point is that the he- 
hx direction in the second order case is determined by 
anisotropics and therefore fixed, whereas it is arbitrary 
in case of the first order transition. In view of MnSi, wc 
propose that this fixed point is the quantum critical end 
point 19] of the observed {p,T) = (12kbar, 12K) tricriti- 
cal point |2i] and might be reachable if one varies another 
thermodynamic variable like chemical composition. 

The equation of state was obtained under the assump- 
tion of a single ordering direction of qp. However, follow- 
ing Ref . [lol one can study more complex choices, like 
a sum over helix configuration with varying direction of 
qo. In case of large number of distinct directions (i.e. an 
amorphous superposition of helix domains) the equation 
of state changes and we obtain, just like in the classical 
case, that these amorphous configurations occur first, at 
least in form of metastable solutions. This is consistent 
with the recent results of Ref.ji^l who showed within a 
replica mean field approach to a similar model that the 
low temperature state might not be perfectly ordered but 
rather be characterized by a distribution of defects of the 
perfectly ordered state with glassy properties. Thus in 
the regime of the I'^^-order transition the system becomes 
sensitive with respect to the smallest amounts of random 
field disorder and forms a self-generated glass. 

If the low energy physics is dominated by the quan- 
tum tricritical point for 

we can also study the effects on fermions scattered by 
magnetic rotons. Since the phase volume of magnetic 
rotons is large, the non-Fermi-liquid properties are even 
more dramatic than close to a 2"'^-order quantum phase 
transitions 0, |2^. Using a local contact coupling with 
coupling constant g, we find the frequency dependence 
of the self-energy of fermions jl^ 

kj, {aEp) ' \EfJ 

(up to a numerical factor), whereas the momentum de- 
pendence of the self energy and the corrections to the 
fermion-roton vertex are nonsingular. The absence of 
singular vertices and the existence of a small coupling 
constant are the key reasons why the the non-Fermi liq- 
uid nature of the fermions does not cause feedback onto 
the collective roton mode, see also Ref. |2^ . Note, Eq|^ 
was previously obtained in Refi for z = 2. Using the 
variational approach to the Boltzmann equation or alter- 
natively a quantum transport theory gives for the resis- 
tivity the result p oc T^' , where 

poc J J d^kd^k'Tk^ (1 - cos^k.k') (10) 

is determined by the angle ^k,k' between k and k' and 
the scattering matrix Tk^k'i determined by Imxq.ujisee 
Ref.jlTjl. Away from the quantum tricritical point, this 



behavior is only expected at higher temperatures, kT ^ 
A. 

Magnetic rotons - magnetic excitations with a min- 
ima on a hyper-sphere - were recently directly observed 
in a weak itinerant helical ferromagnet MnSi by elastic 
neutron scattering^]. In agreement between theory and 
experiment the magnetic phase transition is weakly l'''- 
order type in the range of pressures close to critical pres- 
sure and 2-order type away from critical pressure. The re- 
sistivity in the disordered high-temperature phase shows 
\/r-behavior [2^ . which is consistent with our result us- 
ing z — 2 for the dynamic exponent due to conventional 
Landau damping. 

In summary, we have shown that helicoidal magnets 
with weak anisotropy of the helix direction display a rich 
spectrum of interesting properties. We showed that the 
system is governed by a quantum tricritical point with 
accompanied non- Fermi liquid behavior of the electrons. 
We propose that the (dis)ordered state is likely charac- 
terized by an amorphous superposition of domains with 
different helix directions. 
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